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ABSTRACT 


The primary objective of this program was to demonstrate the feasibility 
of replacing existing state-of-the-art analog gyro compensation loops with 
digital computations. This objective was realized during the course of the 
program. 

A breadboard design was established in which one axis of a Teledyne tuned- 
gimbal TDF gyro was caged digitally while the other was caged using con- 
ventional analog electronics. The digital loop was designed analytically to -v* 
closely resemble the analog loop in performance.' The breadboard was subjected* 
to various static and dynamic tests in order to establish the relative stability 
characteristics and frequency responses of the digital and analog loops. Several 
variations of the digital loop configuration were evaluated. The results were very 
favorable - it appears that digital caging is indeed a practical approach. 
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I. INTRODUCTION 


Background 

In January 1974 Teledyne Systems Company completed under contract to 
NASA Langley Research Center an “Investigation of the Application of Two- 
Degree -of- Freedom Dry Tuned- Gimbal Gyroscopes to Strapdown Navigation 
Systems"*. During the course of that study an entirely new approach to the 
control and compensation of the TDF gyro, as well as accelerometers, was 
suggested and studied. This technique was based upon the use of "all-digital" 
sensor compensation loops and control functions. 

An analysis of the hardware content and costs of present day inertial 
navigation systems showed that a significant portion, on the order of 30%, 
of the system, consists of analog and computer interface electronics. The 
advent of newly available, versatile, digital microcomputers and high-speed 
analog- to-digital and digitai-to-analog converters now makes it practical to 
replace essentially all of the analog electronics by functionally equivalent 
digital computations. 

Figure 1 shows a functional block diagram of a conventional strapdown 
inertial navigation system mechanization. Both the sensor compensation 
functions and the spin motor and pickoff excitations are generated by analog 
electronics which operate independently of the digital processor. The sensor 
torquing currents are fed through precision resistances to develop voltages 
which are proportional to angular rates and accelerations. These voltages 
are then converted into digital numbers for use in the navigation and attitude 
equations which must be solved by the computer. 

^ shows the corresponding block diagram for a strapdown system 
employing digital sensor control and compensation. Here all of the compen- 
sation loop servo functions are mechanized as digital computations in the proc- 
essor. The basic inputs to the A/D converter are the sensor pickoff signals 
rather than the restoring signals. The torquing signals are computed as digital 
numbers, converted to analog signals and, after power amplification, used to 
restore the sensors. These same digital torquing quantities may be used di- 
rectly in the strapdown attitude and navigation computations as well. (In an 
optional configuration the actual torquing signals are also passed through the 
A/D converter for use in the attitude and navigation equations. ) 


*The final report of this study was published as NASA CR-132419 
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Figure 1* Functional Block Diagram of Conventional Strapdown Mechanization 



Figure 2, Functional Block Diagram of '’All-Digital 11 Strapdown Mechanization 
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What is gained by such an “ail-digital" approach is that the presently 
cumbersome and necessarily specialized analog and digital interface elec- 
tronics can be virtually eliminated, being replaced by digital computations 
(software) performed in the computer hardware which is, in any case, 
required for the attitude and navigation computations. This results in sub- 
stantial savings in cost, size, weight, and power consumption while simul- 
taneously increasing the reliability of the system by reducing the component 
count. 


Objective 

The objective of this program was to determine the feasibility of replacing 
present state-of-the-art analog circuits in strapdown tuned-gimbal gyro com- 
pensation circuits with their digital counterparts. This was accomplished by 
designing appropriate compensation loops for the dry tuned TDF gyro, select- 
ing appropriate data conversion and processing techniques and algorithms, and, 
using existing laboratory equipment, breadboarding the design for laboratory 
evaluation. 

The principal area of engineering design involved in the program was the 
determination of the specific software requirements for closing the instrument 
loops and performing the required compensation. Specific requirements rel- 
ative to processing rate, word length, computer time and memory utilization 
were established in order that corresponding requirements for a full comple- 
ment of strapdown sensors may be readily extrapolated. Additionally, appro- 
priate analog-to-digital and digital-to-analog conversion designs were estab- 
lished and techniques selected which were capable of meeting these requirements. 

An available dry tuned-gimbal two-degree-of-freedom gyroscope was util- 
ized together with an electronic circuit breadboard and existing laboratory 
hardware to implement the digital control loop. Selection;of hardware and 
techniques was based on the anticipated availability of proven hardware in the 
late 1970's. The breadboard design was evaluated using standard test and 
evaluation procedures. 

This final report summarizes the test configuration and the program 
results which were obtained. 
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IX. ANALYSIS 


The breadboard design which was established for this study, consisted of 
a dry-tuned two-degree-of-freedom (TDF) gyroscope with one axis being 
caged in a conventional manner using simplified existing "state-of-the-art" 
analog compensation. The second axis was caged using a digital compensation 
loop designed to emulate the performance of the analog compensation. Figure 
3 shows a simplified block diagram of the breadboard design. An analysis of 
the two caging techniques employed in the design is presented in the following 
paragraphs. 


Analog Compensation Design 

The general Teledyne gyro caging loop configuration for conventional 
analog compensation is shown in the simplified block diagram of Figure 4, 
Observe that both direct and cross axis compensation is employed in this 
design in order to increase the loop bandwidth and reduce the rotor hangoff 
during acceleration inputs. In this block diagram and 0y represent fche 
gyro pickoff angles, T x and Ty the gyro torques, F x and P y the pre-. 
cessional torques resulting from rate inputs u> x and w y, and M x and My the 
rebalance torques provided by the compensation loops. 

A more detailed block diagram of the compensation portions of the loop 
for the current Teledyne design is shown in Figure 5. This block diagram 
shows the actual transfer functions which are implemented in the analog 
electronics as well as the pickoff and torquer transfer characteristics. A 
derivation and analysis of this compensation design is included in NASA 
CR-132419. The pole-zero locations for the compensation electronics using 
this design are shown in Figure 6. 

Concurrently with the digital loop demonstration program Teledyne has 
been involved in an extensive redesign effort on its conventional analog 
compensation electronics. This effort has as its goals improved loop per - 
formance and simpler loop design. Although work is continuing in this area, 
it has been found that excellent performance can be obtained with relatively 
simple analog transfer functions. The pole -zero configuration for one such 
simple loop is shown in Figure 7. This configuration provides a baseline for 
much of the digital loop design work which was performed during the program. 
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Figure 3. Block Diagram of the System Used in the Digital Gyro-Caging Study 
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Digital Compensation Analysis 


Figure 8 shows a simplified block diagram of the digital gyro caging 
loop. This figure should be compared with Figure 4. One of the primary 
objectives of the program was to design digital filters -which would emulate 
the performance of the existing analog filter design. 


The study of digital filters and digital servos has progressed rapidly in 
the past decade and an extensive literature now exists with regard to these 
subjects. Rabiner and Rader*, for example, provides both theoretical and 
practical coverage of digital filtering as well as extensive bibliographies. 

The large variety of digital filters and the subtleties (e. g. , the aliasing effect) 
which arise in their implementation, however, serve to maintain digital 
filtering as more of an art than a science. 

Some "cookbook" approaches exist for obtaining digital equivalents to 
analog filters. Some of these have been analyzed, aided where required by 
simulation, in order to determine their applicability to the digital compen- 
sation loop design. In selecting a practical approach, particular attention was 
paid to the impact on the computer with regard to memory and time utiliza- 
tion. As an example, it was desirable to synthesize digital filters whose 
coefficients are integer powers of two, since this allows the replacement of 
multiply instructions with less time consuming shift instructions. 

As one simple example of using the cookbook approaches, consider the 
bilinear transformation approach using the transformation 


s 


1 -z 


1 +z 


-1 

IT 


The transfer function 


H( s) = 


. 79N 


+ 1 


( — + 1) 2 
3N } 


appearing in Figure 5 is seen to have a zero atu^ = . 79N = 496 rad/sec 
and a double pole at 05^ = 3N = 1884 rad/ sec. In order to synthesize a 
digital filter with these characteristics we first derive an analog filter with 
a zero at 


*Rabiner, L. R. and C.M, Rader, Digital Signal Processing, IEEE Press, 
New York, 1972. 
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4C T 

’ 1 
“ i = Un — 


and a double pole at 


* *Z T 

" 2 = tan ~ 


where T is the computer iteration interval. Assuming, for example, an 
iteration rate of 200 Hz, r is 5 milliseconds and 

(496)(. 005) 


to = tan " — 1 — ~ T ■' = 2.91 
1 ^ 


' , (1884)(. 005) _ 

to _ - tan z = 418, 59 


so that the ''dummy” analog transfer function is 

s 


H’(s) = 


2.91 


■ + 1 


S :+l i2 


1 


^418. 59 

Finally, the digital transfer function is found by introducing the transformation 

! -1 
I - z 

s — ► — so that 


1+z 


-1 


H(z) = 


1 d-z" 1 ) 

2.91 (l+z* 1 ) 


+1 


1 (1-z -1 ) 

418. 59 (1+z" 1 ) 


+1 


r 1 n *- 1 \ 4. (U 'M 

[l+z- 1 ] 

L2. 91 (1_z ) + (1+z >J 

•4 

1 

1— * 
J 

t 

►-* 

2 

|^418. 59 1_Z > + ' 
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■ — S— 

418.59 (418.59) 


H' 


(418.59)' 


) 


-1 


z 



2 

(418. 59) 


•f 


1 

(418.59) 



-2 

z 


The bracketed groups in this expression are fixed constants. Thus an 
exact mechnaization of this filter requires 6 multiplies and 4 adds per iteration 
of the computer. Observe, however, that one of these multiplications is by 
the number two. Thus for this operation a shift may replace a multiply ins true 
tion resulting in a substantial saving in computer time since a simple shift 
requires only a small fraction of the time needed for execution of a multiply 
in struction. 


Further simplifications may be possible in this transfer function. For 
example, consider approximating the factor (1 + 1/2.91) = 1.34 by 


1.375 = 1 4*4" -Hr 
4 8 


Then the multiplication (of the current input) by this factor may be accom- 
plished by a double shift, a single shift, and two additions, again effect- 
ing a substantial saving in computer time. Such approximations must be 
carefully analyzed, however. Working "backward" through the preceding 
derivation it is seen that this approximation has the effect of shifting the 
zero of the transfer function from 40^ = 496 rad/ sec to 

2 -1 

10 = — tan (2. 67) 

1 T 


2 

.005 


(1.21) = 485 


rad 

sec 


In this case the effect on the filter response is probably negligible. 
This, is not the case for all approximations which appear on the surface to 
be reasonable, however. 
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Initial efforts in designing the digital compensation filters concentrated 
on emulating existing "old design n analog filters as used in previous strapdown 
systems. The transfer functions of these analog filters are shown in the block 
diagram of Figure 4 and their pole-zero locations in Figure 5. 


Consider first the cross axis transfer function relating My to i. e. 
the y-axis torque r output to the x-axis pickoff angle. The analog transfer 
function is 


M y {s) 
0^4 s) 


'- 01 K sa K po K T 



( 16N + 15 ( 3N +1) 



Since the pickoff and torquer transfer functions will exist in the digital 
loop in any case, they do not have to be implemented digitaly. The gain 
factor (1.01 K ga ) may also be neglected in the digital design as the gain may 
be provided by the power amplifier. Thus the digital transfer function to be 
mechanized should be equivalent to 


H(s) 



Using the bilinear transformation technique and assuming now an 800 Hz 
update rate, the dummy analog breakpoints are computed 


to' 


_ (2 ft) (35) _ _. Q ,_ 

tan (2) (800) # 13832 


to 1 . 


- ■ n (2 ID (79) _ 

(2) (800) " • 32058 


40'. 


= tan (2ff) (100 > g 41421 

an (2) (800) ,4A4 ' 2A 


to' 
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Making the substitution 


1 - z 


-1 


1 + z 


-1 


the digital transfer function is found (after much arithmetic manipulation) 
to be 


H( z ) = 57. 8787 


1 -1 . 68566z ^ +. 5Q190z 
1 +. 65685 z _1 ~.62742z' 2 


41 . 36534 z ~ 3 - 1 , 43278 z ~ 4 + 3 . 8945 z 
- . 74517 z ~ 3 -. 25483 Z -1 -. 02944 z " 5 


In order to obtain a rough estimate of computer timing requirements for 
implementation of the digital loop, the numerator of H{z) was considered. The 
numerator was first approximated in order to eliminate the need for multiplica- 
tions. Thus 

1 . 68566 ^ 1 4 - 2 -1 4 - 2~ 3 + 2 -4 = 1 . 6875 
. 50190 = 2 "* 1 = . 500 


1. 36534 = 1 4 2.“ 2 .+ 2 ~ 3 = 1 . 375 
1 . 43278 ^ 1 + 2 -2 + 2~ 3 -f 2" 4 = 1.4375 
. 38945 = 2 ** 2 + 2“ 3 + 2 -6 = . 390625 


so that the numerator is implemented as 


9 (k) - (l+ 2 ' 1 + 2 ' 3 +E' 4 ) 9 (k— 1 ) + (2 _1 )e (k- 2 ) + (l+ 2 _2 + 2 ' 3 ) 9 (k- 3 ) 

x X X X 

- { 1 +2 "* 2 +2'" 3 +2 ~ 4 ) e (k-4) + (2“ 2 +2“ 3 +2 -6 ) 8 (k-5) 

x x 

(Working "backward” as before, these approximations resulted in analog 
breakpoints of .519N, .749N and . 711N 706N, a reasonable agreement with 

the original transfer function. ) Using the TDY-52B/IMP-16 instructions, a 
simple program was written to implement the equation above. The instruction 
breakdown was as follows 
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In struction 

Number 

Time Per 
Instruction 

Time 
(p sec) 

Load 

6 

7.7 

46.20 

Store 

5 

9. 1 

45.50 

Reg. Copy 

1 

8.75 

8. 75 

Reg. Add 

10 

4. 55 

45.50 

Add/ Sub 

5 

7.7 

38. 50 

Shift Rt. 1 

3 

10. 1 

30.30 

Shift Rt. 2 

5 

14. 3 

71.50 

Shif t Rt . 3 

1 

17. 5 

17.50 




303 . 7 5 psec 


Assuming that each transfer function numerator and denominator is of 
roughly the same complexity, the total time required to impement the four 
transfer functions is 


(303.75) x Z x4 = 2,43 msec 


which is approximately twice the available time for an 800 Hz update rate. 

Since this estimate does not include I/O processing, scaling, executive, con- 
trol, etc. it was clear that much simpler processing would be required if the 
digital compensation loop were to be successfully implemented with the ■ 
TDY-52B computer. 

Fortunately, as described in Section II, efforts were underway concurrent 
to the digital loop program to simplify and improve analog compensation 
techniques. This work resulted in the designs described in the preceeding 
section, and showed that very simple analog transfer functions could be used 
to obtain the desired servo loop response characteristics. Furthermore, it 
was decided to mechanize one axis of the gyro in an analog loop in order to 
more accurately compare the analog loop and digital loop responses, thus 
reducing the digital computations which were required.. These two factors 
resulted in a digital loop design which could be implemented with a TDY-52B 
even though the update rate was increased to 1200 Hz, 

The ’’first cut” at implementing the new transfer functions started with 
the analog functions whose pole-zero locations are shown in Figure 7. 


F^s) 

F 2 (s> 


s + N 
s + 3N 

s 4- N 
s (s + 3N) 
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where F^(s) is the direct axis transfer function and F 2 (s) is the cross axis 
transfer function. Using a 1200 Hz update, the dummy breakpoint frequen- 
cies were computed to be 


w 'i= tan lf?iS = - 26796 


W ' 2 = 




so that 


i z ,~. 1 — j. 26795 

F l(z )= { z+ i! • = - .73205 = . 63 3 98(1 _. 57735z -l ) 


(z-1) 

(z-1) 

(z-1) 

(z + 1) 


+ 1.0000 


+ .26795 


2 z 


r! “ ■ {58 m * ‘-H 


= .63398 


(1 + . 42265z~ 1 - . 57735z~ 2 ; 
(1-z *) 


Since the constant gain term can be readily implemented in the analog 
portion of the loop, it may be neglected. The transfer functions were further 
simplified to 

G^z) = (1 - . 5z 1 ) 


G 2 (z) = 


(1 + . 5z~ ^ .. . 5z" 2 ) 
(l-z -1 ) 


for easier digital implementation. This has the effect of changing 
W'j to .333 and (the analog loop zero) to 1.23N instead of N. 


The comparison of analog and digital frequency responses is shown 
in Figures 9a and 9b for this configuration. It is seen that the breakpoints 
and low — frequency asymptotics are identical, with variations occurring at 
intermediate points. The high frequency performance of the digital loop is, 
of course, different from the analog response due to the "fold-over" effect. 
For this reason, the filter update rate must be sufficiently high that the 
deleterious effects of this phenomenon are countered by attenuation in the 
analog portions of the loop. (Test results indicate that 1200 Hz is the min- 
imum update frequency required to obtain the desired response and that 
higher frequencies are desirable. ) 
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Figure 9a. Direct Axis Frequency Response 
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• DIGITAL LOOP RESPONSE 
O ANALOG LOOP RESPONSE 


Several variations of parameters were tested within the constraints 
of the basic filter configurations described by G;[{z) and G£ (z) above. 

It was found that the transfer functions 


G x (z) 
G 2 (z) 


- n 15 - 1 ! 

~ (1 "l6 > 

1 + —z~^ 

.16 2 


1 -z 


-1 


15 

16 


-2 

z 


provided reasonably optimal performance for the digital loop*. Most 
of the results described in the report reflect this configuration. Thus the 
equations which were mechanized in the computer were 

M (k) = M (k-1) + e (k) + (2 ” ) e (k-1) - (1 -2*" 4 ) 0 (k-2) 

XXX X X 

M (k) = 0 (k) - (l-2“ 4 ) S (k-1) 

y y y 

M(k) = M (k) - M (k) 
x y 

Jt»T_ 

where M(k) is the computed gyro torque at the k l iteration. 

The digital number M(k) is converted by the digital-to- analog 
converter and, after power amplification, used to torque the gyro. In 
the absense of errors, M(k) represents the true torque input and is, therefore 
proportional to the angular rate of the gyro case. This is precisely the 
information which is desired for attitude computation in a strapdown navigator 
Any errors introduced by the D/A converter therefore propagate directly 
into attitude errors. For this reason it is necessary to use an accurate D/A 
converter in the digital loop. 


*T hi s corresponds to moving the zero locations in Figure 6 from 100 Hz 
to approximately 12.7 Hz. 
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The analog-to -digital conversion requirements are much less stringent. 
Due to the closed-loop nature of the digital compensator, errors in the A/D 
converter are of relatively minor importance to the accuracy of rate measure- 
ment (just as pickoff measurement errors are unimportant in a conventional 
analog loop mechanization). 

The test results indicate that a 16 bit computer word length is sufficient 
for mechanization of the digital loop. The D/ A converter should be an accurate 
device with excellent (or, at least, well compens atable) bias and scale factor 
characteristics. The A/D converter can be less accurate without seriously 
compromising torque measurement accuracy. The digital loop was mechanized 
using a 12 bit A/D converter and a 13 bit D/A converter. These resolutions 
are not critical, however, for the reasons outlined above. As long as the 
analog torque applied accurately represents the (quantized) digital torque, no 
net error will result in the computation of attitude due to the' D/A converter 
resolution. 
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III. COMPUTER 


The computer used in the breadboard is an IMP-16P (National Semi- 
conductor), which is electrically and functionally interchangeable with the 
T SC hybrid packaged TDY-52B. Although these two machines are. function- 
ally the same, they are vastly different physically. The CPU of each con- 
sists of four 4 -bit NS (National Semiconductor) GPC/P RALU slices. The 
microprogram of each is contained in a set of two Control Read Only Memories 
(CROMs) that implement the NS IMP-16 computer repertoire. Input/Output 
and other hardware details are also similar. The IMP-16P and TDY-52B 
are thus, basically, the same, machine. 

The most obvious difference between the IMP-16P and TDY-52B is one 
of size. Teledyne has packaged the heart of the IMP-1 6F into a small hybrid 
unit that is only two inches on a side and .2" thick. The IMF-16P computer, 
including a convenient control panel and 4K of 16 -bit RAM memory, occupies 
a box with a frontal dimension of about 10 1/2" x 17" and a depth of about 24". 

The IMP-16P was selected for the digital loop demonstration pri- 
marily because its control panel feature facilitated the operation and 
modification of, the breadboard. The results which were obtained are, 
however, equally applicable to the TDY-52B. 
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Table 1. TDY-52B/IMP-16 Instruction Repertoire (Sheet 1 of 2) 


MNEMONIC 





EXECUTION TIME 

INSTRUCTION NAME 

FUNCTION 


FORMAT 

IN MICROSECONDS 

LOAD AND STORE 





7.0, 9.8 
IF INDIRECT 

LD 

LOAD 

(EA)- (ACr)-ylF INDIRECT ((EA))-(ACr) 

2 

LDB 

LOAD BYTE 

(1/2 EA)-(ACO LESS SIGNIFICANT BYTE) 

5 

16.8 TO 28.0 

LI 

LOAD IMMEDIATE 

' D - (ACr) 


4B 

4.2 

ST 

STORE 

(ACN)- (EA), IF INDIRECT (ACr)-((EA)) 

2 

8.4,11.2 IF INDIRECT 

STB 

STORE BYTE 

(ACO LESS SIGNIFICANT BYTE) — (1/2 EA) 

5 

23.8 TO 32.2 



O-(SEL) 




RXCH 

REGISTER EXCHANGE 

(SR)-* (DR), (DR)-* (SR) 


6 

11.2 

RCPY 

REGISTER COPY 

(SR) -*(DR) 


6 

8.4 

LOGICAL 






RXOR 

REGISTER EXCLUSIVE OR 

(SR) ©(DR) * (DR) 


6 

8.4 

RAND 

REGISTER AND 

(SR) “AND" (DR)-* (DR) 


6 • 

8.4 

AND 

AND 

(ROD "AND" (EA)--(ROl) 


3 

7.0. . 

OR 

OR 

(R01) "OR" (EA)-* (R01) 


3 

7.0 

ARITHMETIC 






RADD 

REGISTER ADD 

(SR) + (DR)-*(DR) OV, CY 


6 

4.2 

ADD 

ADD 

(ACr) + (EA)-* (ACr) OV, CY 


2 

■7.0 

SUB 

SUBTRACT 

(ACr) - (EA)- (ACr) OV, CY. 


2 

7.0 

MPY . 

MULTIPLY 

(EA) * (AClHj\AC0, lACl)] L O-(SEL) 

5 

148.4 TO 170.8 

DIV 

DIVIDE 

( (ACO), (ACl)j-(EA)— (ACO) QUOTIEN7 
0 - (SEL) OV, L Nacd REMAINDER 

5 

177.8 TO 222.6 

DADD 

DOUBLE PRECISION ADD 

((ACO), (ACI)| + j(£A), (EA+])1-HAC0), (AC 1)1 
0— (SEL) OV, CY 1 11 1 

5 

16.8 

DSUB 

DOUBLE PRECISION SUBTRACT 

|(AC0), fACl)j- j(EA), (EA + 1)U((AC0), (ACI)i 

5 

16.8 

CAI 


0 — * (SEL) OV, CY 


4B 

4.2 

COMPLEMENT AND ADD 

~ (ACr) + D- (ACr) . 




SHIFT 

IMMEDIATE 




ROL 

ROTATE LEFT 

2 (ACr) - (ACr) 

IF SEL - 0, (BIT 15)— (BIT 0} 

IF SEL = 1 , (BIT 15) — (L) (L)-(BIT 0) 

| D TIMES 

4B 

5.6*4.20 


ROR 

ROTATE RIGHT 

1/2 (ACr) — (ACr) 

IF SEL ~ 0, (BIT 0) — (BIT 15) 

IF SEL = f, (BIT 0) — (L) , (L)-(BIT 15) 

D TIMES 

4B 

5. 6 + 4. 2D 

SHL 

SHIFT LEFT 

2 (ACr) - (ACr) 

0— (BIT 0) 

IF 5EL = 1 , (BIT 15)- (L) 

D TIMES 

4B 

5.6 +4. 2D 

SHR 

SHIFT RIGHT 

1/2 (ACr)— (ACr 1 
IF SEL = 0, 0-(BIT 15) 

IF. SEL = 1, (L) -(BIT 15), OWL) 

D TIMES 

4B 

5.6 + 4. 2D 

SINGLE BIT 






SETST 

SET STATUS BIT 

1 -{STATUS FUG N) 


9 

18.2 TO 44.8 

CLR5T 

CLEAR STATUS BIT 

0- (STATUS FUG N) 


9 

18.2 TO 44.8 

SETBiT 

SET BIT 

1— (ACO BIT N) 


9 

18.2 TO 44.8 

CLRBIT 

CLEAR BIT 

0— (ACO BIT N) 


9 

18.2 TO 44.8 

CMPBIT 

COMPLEMENT BIT . 

(ACO BIT N) — (ACO BIT N) 


9 

18.2 TO 44.8 
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Table 1. TDY-52B /IMP- 1 6 Instruction Repertoire (Sheet 2 of 2) 


MNEMONIC 

INSTRUCTION NAME 

FUNCTION 

FORMAT 

EXECUTION TIME 
IN MICROSECONDS 

JUMP 





JMP 

JUMP 

EA-(PC) , IF INDIRECT (E A)— - (ACr) 

4A 

4. 2,7.0 IF INDIRECT 

JMPP 

JUMP THROUGH POINTER 

(100 )6 +N)- (PC) 

9 

9.8 

JINT 

JUMP INDIRECT TO LEVEL 0 

(PC) — (STIC), O-(IEF) 

9 

9.8 



(120 |6 +N)- PC 



BOC 

BRANCH ON CONDITION 

IF CONDITION CC IS TRUE, 
(PC) + D - (PC) 

1 

5. 6, 7.0 IF BRANCH 

J5R 

JUMP TO SUBROUTINE 

(PC)-(STK) 

EA — (PC) , IF INDIRECT 'EA) — (PC) 

4A 

5. 6,8. 4 IF INDIRECT 

JSRI 

JUMP TO SUBROUTINE IMPLIED 

(PC) - (STK) 

8 

5.6 



FF80, + C-(PC) 



• JSRP 

JUMP TO SUBROUTINE THROUGH 

f PC)—'STK) 

8 

11.2 


POINTER 

(100 t6 +C)-(PC) 



RTS 

RETURN FROM SUBROUTINE 

(STK) +C-(PC) 

8 

5.6 

RTI 

RETURN FROM INTERRUPT 

(STK) + C — (PC) 

8 

7.0 



l-(IEF) 



SKjP 





AISZ 

ADD IMMEDIATE AND 
SKIP IF ZERO 

(ACr) + D — (ACr) OV, CY 
IF (ACr) = 0, (PC) + 1— -(PC) 

4B 

5. 6,7.0 IF SKIP 

SKAZ 

SKIP IF "AND" IS ZERO 

IF (ROI) "AND" (EA) - 0, (PC) i 1— (PC) 

3 

8.4, 9.8 IF SKIP 

ISZ 

INCREMENT AND SKIP IF ZERO 

(EA) + 1— (EA) 

IF (EA) = 0, (PC) + 1 — (PC) 

4B 

9.8, 11.2 IF 5KIP 

OSZ 

DECREMENT AND SKIP IF ZERO 

(EA) - 1 — (EA) 

IF (EA) - 0, (PC) + 1 - (PC) 

4B 

11.2,12.6 IF SKIP 

SKG 

SKIP IF GREATER THAN 

IF (ACr)»(EA), (PC) + l-(PC) 

2 

11.2 TO 14.0 

SKNE 

SKIP IF NOT EQUAL 

IF (ACr) 4 (EA), (PC) + 1 - (PC) 

2 

8.4 

skstf 

SKIP IF STATUS FLAG TRUE 

IF (STATUS FLAG N) = 1, (PC) + 1— (PC) 
0 (SEL) 

9 

18.2 TO 44.8 

SKBIT 

SKIP IF BIT TRUE 

IF (ACO BIT N) = 1, (PC) + 1 — (PC) 
0 (5EL) 

9 

18.2 TO 44.8 

STACK 





PUSHF 

PUSH STATUS FLAGS 
ONTO STACK 

(SF) —(STK) 

8 

5,6 

PU.LLF 

PULL STATUS FLAGS FROM 
STACK INTO FLAG REGISTER 

(STK)- (ACr) 

8 

7.0 

PUSH 

PUSH ONTO STACK 

(ACr) -(STK) 

4B 

4.2 

PULL 

. PULL FROM STACK 

(STK)- (ACr) 

4B 

4.2 

XCHRS 

EXCHANGE REGISTER 
AND STACK 

(ACr) —(STK) 
(STK)- (ACr) 

4B 

7.0 

input/output 





RIN 

REGISTER INPUT 

(AC3) + C-(IO ADDR) 
(IO DATA) -(ACO) 

a 

9.8 

ROUT 

REGISTER OUTPUT 

(AC3) + C - (IOADDR) 
(ACO)- (10 DATA! 

8 

9.8 

SFLG 

SET FLAG 

C— OOADDR) , 1- (CONTROL FLAG FC) 

7 

5.6 

PFLG 

PULSE FLAG 

C— ((OADDR), I— (CONTROL FLAG FC) 

7 

5.6 

0 ISCAN 

INTERRUPT SCAN 

1/2 (ACl)-(ACl) UNTIL 1 SHIFTED OUT 
<AC2) + NUMBER OF SHIFTS-(AC2) 

9 

8.4 TO 100.8 

HALT 

HALT 

PROCESSOR HALTS 

8 

' 
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Notation Used in Instruction. Descriptions 



Denotes a specific working register (ACO, AC1, AC2, or AC3), where r is the number of the accumula- 
tor referenced in the instruction. 


Denotes the address register used for addressing memory or peripheral devices. 


Denotes the 4-bit condition code value for conditional branch instructions. 







Denotes the 7-bit control-field value for flag, input/output, and miscellaneous instructions. 


Indicates that the Carry flag is set if there is a carry due to the instruction (either an addition or a 
subtraction). 


Stands for displacement value and it represents an operand in a nonmemory reference instruction or an 
address field in a memory reference instruction. It is an 8-bit, signed twos-complement number except 
when base page is referenced; in the latter case, it is unsigned, 


Denotes the number of a destination working register that is specified in the instruction-word field 
The working register is limited to one of four: ACO, ACl, AC2, or AC 3. 


Denotes the effective address specified by the instruction directly, indirectly, or by indexing. The 
contents of the effective address are used during execution of an instruction. See table 3-1 . 


Denotes the number of the referenced flag (see table 3-20 under 3.6.10, Input/Output, Halt, and 
Flag Instructions). 


Denotes the Interrupt Enable control flag. 


IOREG Denotes an input/output register in a peripheral device. 


Denotes 1 -bit link (L) flag. 


Indicates that the overflow flag is set if there is an overflow due to the instruction (either an addition 
or a subtraction). 


Denotes the program counter. During address formation, it is incremented by 1 to contain an address 
1 greater than that of the instruction being executed. 


Denotes the number of a working register that is specified in the instruction-word field. The working 
register is limited to one of four: ACO, ACl, AC2, or AC3. 

Denotes the Select control flag. It is used to select the carry of overflow for output on the carry and 
overflow (CYOV) line of the CPU, and to include the link bit (L) in shift operations. 


Denotes the number of a source working register that is specified in the instruction-word field. The 
working register is limited to one of four: ACO, ACl, AC2, or AC3. 


When not zero, this value designates the number of the register to be used in the indexed and relative 
memory-addressing modes. 
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Notation Used in Instruction Descriptions (Continued) 


Notation 

Meaning 

( ) 

Denotes the contents of the item within the parentheses. (ACr) is read as “the contents of ACr 
(E A) is read as “the contents of EA “ 

1 1 

Denotes “the result of.” 


Indicates the logical complement (ones complement) of the value on the right-hand side of 


Means “replaces.” 

+- 

Means “is replaced by.” 


Appearing in the oparand field of an instruction, denotes indirect addressing. 

A 

Denotes an AND operation. 

V 

Denotes an OR operation. 

V 

Denotes an exclusive OR operation. 
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IV. SOFTWARE 


Introduction 

The intent of this study was not to demonstrate the particular processor 
which was used but rather to demonstrate how the program was mechanized 
for this particular computer in such a manner that the software could be 
mechanized for any given computer with equivalent characteristics. The 
reader should have some knowledge, however, of the computer in which the 
equation set was mechanized in order to fully understand the work performed. 
An example of how the flow charts reflect the instruction repertoire of the 
processor is in the handling of negative quantities. The TDY-52B/IMP- 16 
has a logical shift only, so that when a multiplication by 1 / 1 6 is to be per- 
formed it is necessary to handle negative quantities in a different manner than 
positive quantities where the multiplication can be performed by a shift of 4 
bit positions. 

In order to give a better understanding of the digital computer used in 
this study a summary of instructions was included in Section III. Should the 
reader have further interest in this processor including the interrupt system, 
push-pull stack operation, logic mechanization, etc, , complete descriptive 
manuals can be obtained from Teledyne Systems Company. 

Power On Processing 

Upon recognition of Power being applied to the unit, coding is performed 
to disregard the initial multiplexer interrupt. This is done to allow the pro- 
gram to become synchronized with the timing base of the multiplexer so that, 
after the initial interrupt the timing base of approximately 104 microseconds 
between interrupts is established for the remaining processing. The next task 
is to. initialize a pointer within the processor so that whenever data is gener- 
ated for storage the starting point for this storage is established. 

The interrupts are then enabled so that the processor is ready to recog- 
nize interrupts from the multiplexer. Finally a transfer is made to the Main 
Processing Loop. As the processor used in this study turns on with Random 
Access Memory (RAM) set to an all one’s condition, the initial iterations of 
the Main Processing do not have the proper data for use as a priori iteration 
data so that the data generated for the first two iterations of the Main Process- 
ing Loop is not precise. 


27 



POWER ON INITIALIZATION 



798824 


Main Programming Loop 

The Main Programming Loop is where the majority of calculations are 
performed. This routine is coded to run continuously, i. e. from Power On 
the calculation process is repeatedly performed with the only interruptions 
being the X peak and Y peak Interrupts where data from the analog to digital 
converter is input. Once the data has been input within either of these 
interrupt routines, program control is transfered back to the Main Program- 
ing Loop at the point the program was interrupted. 


Essentially the processing within this software module consists of solving 
the equations below: 

M x (K) - M X (K-1) + U X (K) + 1/16 U x (K-l) - 15/16 U x (K-2) 

M y (K) = Uy (K) - 15/16 (K-l) 

M (K) = U X (K) - U y (K) 
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where: 


M x = filtered digital quantity from X axis pickoff 

My = filtered digital quantity from Y axis pickoff 

M = total digital torque output to D/A 

Index (k) designates quantity input or calculated on current program 
iterations. 

Index (k-1) designates quantity input or calculated on previous iteration 
Index (k-2) designates quantity input or calculated two iterations previous 

The other calculations performed in this loop and shown on the flow 
diagrams for the Main Processing Loop is the interrogation of two software 
flags, X peak flag and Y peak flag, to properly correlate the data from the 
storage buffers and to transform the data from its 2’s complement form for 
the sign magnitude digital to analog convertor. 

Following is the flow chart for the Main Processing Loop. (Figure 10.) 
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MAIN PROGRAMMING LOOP 
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Figure 10. Main. Programming Loop 
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MAIN PROGRAMMING LOOP 



1 MPl 


YES 

Y PEAK^^ 

NO 

| FLAG = 7 

FLAG =7 

FLA G = 0 v 

COMPUTE: 

COMPUTE: 

MY <K) = UY (K) - 15/16 UY (K - 7) 
WHERE UY (K) AND - 15/16 UY 
(K - 1) ARE FROM BUFFER A 


MY (K) = UY (fC) - 75/76 UY (K - 7) 
WHERE UY (K) AND UY (K -1) ARE 
FROM BUFFER B 


COMPUTE: 

Mp = MX (K) - MY (K) 


YES 


PRODUCES 

A SIGN 

MAGNITUDE 

QUANTITY 

FROM A2'S 

COMPLEMENT 

QUANTITY 



! 

r MT< 0 


COMPLEMENT 

MT 

•< 

i 

r 


INSERT 
SIGN BIT 






STC 

/v 

3RE: • 
\T 

< 

r 

OUTPUT: 

MJ 



Figure 10. Main Programming Loop (Continued) 
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X Peak Interrupt Servicing 


Upon recognition by the processor that sample data from the X axis pickoff 
is available, the Main Processing Loop is interrupted at the end of the instruc- 
tion it is performing and the program is transferred to the X Peak Interrupt 
Servicing routine. The purpose of this routine is to accept the digitized X 
axis pickoff data for the computation to be performed by the Main Processing 
Loop. In order to keep the data being computed in the Main Programming 
Loop correlated it was necessary to mechanize two buffer areas which are 
merely memory cells set aside for temporary data storage and retrieval. 

At any given point of time one of the buffer areas is set for reading by the 
Main Processing Loop and the other buffer area is set for writing or storage 
by the X Peak Interrupt Servicing routine. The purpose of these buffers is 
to ensure that the data being computed is indeed from successive iterations. 

A memory cell is used as a flag (X Peak Flag in flow charts) to indicate how 
the buffers are currently being utilized. This s'ame flag is again interrogated 
in the Main Processing Loop for the same purpose. A graphic illustration 
of the storage areas referred to as ’swinging buffers* will perhaps clarify 
this coding concept. At the time the X Peak Interrupt is entered by the 
program assume that the two buffers are as shown below 


Buffer A 


Buffer B 


and further assume that the flag (X Peak Flag) indicates that Buffer A is the 
area that should be used to store data. The cell from Buffer B containing 
UX(K-l) must first be moved to Buffer A as it now becomes two iterations 
old, i. e. U x (K-2). In the same manner the Buffer B cell containing Ujr(K) is 
moved to Buffer A and is now one iteration old, i. e. U x (K-l). The digitized • 

X axis pickoff data input during the current iteration is then stored in Buffer A 
as Ux(K). Buffer A then contains data from three successive inputs. The 
process is reversed the next iteration such that Buffer B contains current data. 
The last task of the X Peak Interrupt processing is to divide the data from the 
previous iteration by 16 (accomplished by a right shift of 4 places) before 
storing to relieve the Main Processing Loop of this computational burden. 

Once these tasks are performed the program is returned to the Main Process- 
ing Loop, to the point at which the interrupt occurred. 

Following is the flow charts for the X Peak Interrupt Processing. 

(Figure 11.) 


V K > 

U (K-l) 
X 

U x (K-2) 


— 

"x 

D 

U X (K-1) 

U x (K-2> 
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X PEAK INTERRUPT SERVICING FLOWS 



Figure 11, X Peak Interrupt Servicing Flow 
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Y Peak Interrupt Servicing 


Upon recognition by the processor that sample data from the Y axis 
pickoff is available the Main Processing Loop is interrupted at the end of 
the instruction it is performing and the program is transferred to the Y Peak 
Interrupt Servicing routine. The purpose of this routine is to accept the 
digitized Y axis pickoff data for the computation to be performed by the Main 
Processing Loop. This data is correlated in much the same manner as 
previously explained in the X Peak Interrupt Servicing in that a flag (Y Peak 
Flag) is used to determine which data is time correlated and buffers are 
used to determine read and write areas. In this process the last iteration 
Y axis pickoff data Uy(K- 1) is multiplied by 15/16 as a time saving step to 
the Main Processing Loop. Again, as in the X Peak Interrupt Processing, a 
return is made to the point in the Main Processing Loop where processing 
was interrupted. 

Following is the flow charts for the Y Peak Interrupt Processing. 

{Figure 12,) 
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Y PEAK INTERRUPT SERVICING FLOWS 
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STORE: 

-15/16 K VALUE 
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NEGATIVE 


STORE: 

-15/16 K VALUE 
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(K-l) VALUE 
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QUANTITY ^ 
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STORE: 

>15/16 K VALUE 
IN BUFFER B 
(K-l) VALUE 


NO 

NEGATIVE 




STORE; 

-15/16 K VALUE 
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CONVERTED DATA 
INTO K VALUE 
OF BUFFER A 


STORE NEW 
CONVERTED DATA 
INTO K VALUE 
OF BUFFER B 


RETURN TO POINT 
WHERE INTERRUPTED IN 
MAIN PROCESSING LOOP 


-15/16 UY (K-l) = -1 + 1/16 UY (K-l) 


Figure 12, Y Peak Interrupt Servicing Flow 
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V. ANALOG INTERFACE 


The analog interface for the IMF 16 computer consists of an A/D 
converter, a D/A converter, and the control logic for the two converters. 

A block diagram of the interface is shown in Figure 13. 

The A/D converter section, shown in Figure 14, consists of an input 
multiplex switch, a sample and hold, a 12 -bit converter, and a tri-state 
buffer register. The input multiplex selects an input channel which may be 
a peak, zero, or ground input from one of the gyros. This input depends on 
which channel is selected by the computer and logic decode. The input is 
applied to the sample and hold circuit, which will hold the input signal value 
until the A/D performs a conversion. The A/D start command and clock are 
supplied by the control logic. The tri-state buffer holds the A/D conversion 
and outputs it to the CPU until the next input enable allows the data to update. 

The D/A converter accepts a 13-bit sign-magnitude word from the CPU 
and converts it to a -f or - DC output to drive the gyro torquer amplifier. 

This converter consists of a + DC reference, a 12 -bit ladder network, plus 
output buffer amplifiers. The control logic decodes the input data word 
and outputs a sign bit so that the proper polarity is applied to the ladder 
network. The control logic also applies a 12 -bit word to the ladder network 
to determine amplitude. The signal from the ladder network is applied to 
output buffer amplifiers. The control logic enables either a + or - sample 
switch which applies the respective output signal to the hold capacitor located 
on the torquer amplifier on the platform. 

The control logic of the analog interface receives the 5. 72 MHz CPU 
clock and divides it by 10 and 12 successively to provide basic square waves 
of 572 KHz and approximately 48 KHz, The 48 KHz clock is used to generate 
periodically four sequential analog select signals, each spanning an interval 
of 32 clock periods. Two successive analog-to -digital conversions are made 
during an analog select interval, one in the first half -interval and one in the 
second half -interval. The 48 KHz bit times of each half -interval are 
designated Tq through T^. Analog sampling signals are enabled during 
of each half -interval such that an analog peak signal is sampled in the first 
half- interval and an anlaog zero crossing is sampled in the second half- 
interval. Analog to digital conversion is accomplished during T£ and of 
each half cycle. An interrupt request is then sent to the CPU at T 4 . The 
converted 12 -bit analog data may be input to the computer by the command, 
RlN (0070) The analog data is identified by a three-bit code which is 
inputted to the computer by the command, Rj^ (0078) 16 . 
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The CPU outputs a 13 -bit word in sign -magnitude format for conversion 
to analog by the command, JRqUT (0078)j^. The data word is latched in a 1st 
rank buffer from which it is transferred to a 2nd rank buffer at bit time Tq of 
a continuous 8-bit cycle of the 48 KHz square wave. The conversion is 
accomplished during each To Also, sampling of the converted output is 
disabled during Tq. 
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VI. TEST RESULTS 


Test Description 

Figures 15 through 17 illustrate the general test setup. In this test 
configuration, the gyro was mounted on the rate table such that the input 
rates would be about the gyro x axis. 

The x and y gyro pickoff signals were fed to both the digital servo and 
the analog servo. The digital servo controlled only the y torquer, and the 
analog servo controlled only the x torquer. 

The rate table was driven by the oscillator portion of the wave analyzer 
with a constant voltage applied to the rate table. 

The wave analyzer was used to measure the torquer signals and 
pickoff signals. The x axis (direct) pickoff signal was also measured in the 
open loop condition. The open loop and closed loop measurements were 
used to calculate the gains that are plotted in Figures 18 and 19. 



T97308 


Figure 15. Block Diagram, Test Set Up 
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Figure 16. Test Set-Up for Digital Compensation Loop Feasibility Study 
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Figure 17. Test Set-up of Converters and Interface for Digital 

Compensation Loop System 






OUTPUT SIGNAL (dB) 



ZERO dB = OUTPUT AT DC RATE INPUT 


DATA TAKEN 8.2.74 TMW T97307 

Figure 18. Output Signal (Torquer Voltage) vs Frequency 

Digital -Analog Loop — 100 Hz Bandwidth 









The rate table drive frequency was varied from 10 Hz to 1 KHz with 3 
to 5 measurements taken per decade depending on how rapidly the data was 
changing with frequency. 

The integral gain time constant in the digital servo was modified 
(software modification) and the test repeated. There were no modifications 
to the analog servo. 


Discussion of Test Results 

The test setup described above provided a means for one-to-one 
comparison of the performance of the digital and analog caging loops. In 
order to obtain such a comparison, however, it was necessary to match the 
gains of the two loops since the digital computations, as described in Section 2, 
assumed that the gain would be adjusted on the analog level, i. e. , in'the power 
amplifier. The criterion used to match the gain (applied torque vs pickoff 
angle) of the digital loop to that of the analog loop was to operate each servo 
6 db below its marginal stability level. Using this criterion it was found that 
the digital loop gain was approximately 10 db below that of the analog loop. 

Frequency response data was taken for both the digital and analog loops. 
Figure 20 shows the output signal (i. e. , the torquer voltage) response versus 
frequency. Figure 21 shows the error signal (i. e. , pickoff angle) versus 
frequency. These figures are actually families of responses showing varia- 
tions resulting from adjusting the "integral gain time constant" in the digital 
loop. This time constant is directly related to the location of the zero of the 
digital transfer functions described in Section 2. The range or variation 
corresponds to varying this zero location between 123 Hz and 12. 7 Hz. 

It is apparent from Figure 18 that the closed loop response of the 
digital loop compares very favorably with that of the analog loop in terms of 
the torquer output versus frequency. This quantity is, of course, of primary 
interest as it is the ultimate measure of gyro rate. 

l^l£>n.re 19 indicates that the gyro pickoff angles are substantially larger 
for the digital loop than for the analog loop for low frequency inputs. This 
is attributable to the lower gain in the digital loop which was established by 
the gain matching criterion described above. 
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The primary reason that the marginal stability gain for the digital loop 
is lower than for the analog loop is the effective phase shift introducted by 
the finite computer iteration rate. Thus as the computational frequency is 
increased it is expected that the digital loop gain established by the gain 
matching procedure would be increased and the pickoff angle versus frequency 
response of the digital loop would more closely approximate the corresponding 
analog response. Unfortunately this conclusion could not readily be verified 
due to computer timing limitations. 

The primary ramifications of the lower loop gain in the digital loop is 
a lowering of the bandwidth of the caging loop which would occur in a gyro with 
both axes mechanized digitally. Subsequent analysis indicates that the effec- 
tive bandwidth of such an all -digital gyro would be approximately 25 Hz for 
the configuration which was tested (i. e. , 1200 Hz digital processing). This 
corresponds closely with the 30 Hz bandwidth used in Teledyne's current 
analog caging loops. This bandwidth can be further increased by increasing 
the computational frequency and can be expected to approach the 100 Hz of 
the new analog design with very high speed processing. 

In summary, then, the test results were in general concert with 
expectations and the feasibility of using digital filtering in place of conven- 
tional analog caging loops for the TDF tuned-gimbal gyroscope was demon- 
strated. 
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VII. summary and conclusions 


The primary objective of this program was to demonstrate the feasibility 
of replacing existing state-of-the-art analog gyro compensation loops with 
digital computations. This objective was realized during the course of the 
program. 

A breadboard design was established in which one axis of a Teledyne 
tuned-gimbal TDF gyro was caged digitally while the other was caged using 
conventional analog electronics. The digital loop was designed analytically 
to closely resemble the analog loop in performance. The breadboard was 
subjected to various static and dynamic tests in order to establish the relative 
stability characteristics and frequency responses of the digital and analog 
loops. Several variations of the digital loop configuration were evaluated. 

The results were very favorable - it appears that digital caging is indeed a 
practical approach. 

The primary limitation to the digital compensation loop approach appears 
to be the requirement for high processing rates. The tests which were per- 
formed indicate that the minimum processing rate for acceptable loop per- 
formance is approximately 1200 Hz. Using the general purpose IMP- 16 on 
TDY-52B computers essentially all of the available computer time was re- 
quired to perform the digital caging loop functions for a single gyro axis 
at this rate. (These functions include all timing, control, sampling, etc. , 
as well as the implementation of the actual digital filters. ) For this reason 
it appears desirable to consider a special purpose digital processor for the 
digital caging loop task. Such a processor could be quite simple in concept 
since the only arithmetic functions required to implement the digital filters 
which were derived during the program are adds, subtracts, and shifts. 

The primary advantages of using digital caging loops in place of convex 
tional analog loops are reduction of cost, size, weight, and power, increased 
reliability, and simplified maintenance, resulting from the elimination of 
a substantial amount of hardware from the IMU Electronics. The electronics 
which are eliminated using the digital technique include pickoff demodulators, 
caging electronics, spin .supplies, clock and timing functions, and pickoff 
excitation generation. New functions which must be performed with the digital 
loop are digifcal-to-analog conversion and, optionally, the special purpose 
digital processing. Since both of these functions may be multiplexed to handle 
all sensor axes, the hardware required is minimal. A preliminary design 
analysis has indicated the hardware requirement to perform these two 
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functions is of the order of 5 to 20% of that required for the analog functions 
which they replace depending upon the level of redundancy which is employed 
in the system. The resulting savings in direct acquisition cost for a fail-op/ 
fail-op redundant strapdown navigator is projected to be approximately $4000 
per system in quantities of 150 systems. 

The performance of a system using digital caging loops is expected to be 
virtually identical to that using conventional analog technology. This is 
particularly true if a special purpose processor is employed, since the 
processing rate can be increased beyond the 1200 Hz which was used in the 
demonstration. Higher rates will provide digital loop performance which 
matches the analog loop performance even more closeLy than described in 
this report. 

Other than the limitations imposed by computational speed, the digital 
loop which was breadboarded during the program appears to be suitable for 
use in a practical system design without any significant modifications. The 
general design philosophy, filtering algorithms, computer word length, A/D 
and D/A converter accuracies and resolutions, etc., which were employed 
in the breadboard design are all suitable for use in a production design. 
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APPENDIX A 


PROGRAM LISTING 


This is the assembly language gyro program as assembled on Teledyne 
unspecialized assembler TUMPA, Assembled code was punched out in special 
loading format on mylar teletype tape for loading via teletype into the IMP-16P. 
This program will also run on the TDY-52B. 
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2/OACC, 2/ 

CGC6 

J.PCIF 

2,0SLB,2/ 

C0C7 

A.PCIF 

2,NPY,2/ 

COOb 

• • PC IF 

2,0 IV, 2/ 

COC 9 

J.TRA'i 

3,2/ 

C010 

4. TRAN 

4,2/ 

con 

4 • T ft A N 

3,1/ 

CO 12 

4. TRAM 

4,1/ 

con 

I.TRA r. 

5, 1/ 

CO 14 

*«BOCL 

2, RO, 0/ 

CO 1 5 

4.BMCL 

2, 0 l, 1/ 

con 

i.BOCL 

2, R 2, 2 / 

CQ17 

I.BCCl 

2> R 3, 3/ 

con 

I.BCCL 

2 , acr,c/ 

con 

I.BOCL 

2, AC 1/ I / 

0020 

*.BCCL 

2, AC 2, 2/ 

CO? 1 

i.enci 

2,AC3,3/ 

CC2 2 

j.eccL 

2,nASE,0/ 

C023 

meet 

2, PC, 1 / 

CO 24 

i.bccl 

2,XR2 ,2/ 

CQ25 

i.acu 

2,>t«3,3/ 

CO 2 6 

i.bctl 

2,XMIT,4/ 

0027 

1 , BOC L 

2,LNK,17/ 

CC28 

I.9CCL 

2,0V, 16/ 

CO 2 9 

1 • BCCL 

2,CYi 15/ 

C030 

* < ft n r L 

lfltt.il/ 

C03 1 

I.BCCL 

2,SEL,2/ 

C012 

i.ancL 

2,Rf:ACRi3/ 

con 

A • ti n C L 

?,XH1T,4/ 

C 03 A 

4.H0CL 

2 > JC 14, 16 / 

C035 

I.BDCL 

2, JC 1 , 1 / 

C03b 

t.aoci 

2, JC2, 2/ 

CO 37 

inert 

2, JC3,3/ 
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3.CRC/1/ 



CQ *6 

*.c/ interrupt executive 


0000 1 

C0*7 

* PUSH, *CC/ 


000 C 2 

C 0* 6 

« PUSH* *03/ 


O00C3 

C C *9 

#Ui«3/0/ 


oooc* 

COSO 

*R JN,170/ 


00005 

CO 5 1 

XX,RCPYjACQ,AC3/ 


0000 ft 

C052 

/ JMPiAC3#C/ 



CO 53 

S.CRGMu/ 



cos* 

*.(,/ X PEAK BUFFERS 


000 1 o 

CO 5 5 

XPl,CCT,C/ FLAG1 


ooo n 

C0S6 

XP2,i:Cr#C/ CELL 2 


00012 

C057 

XPi#0CT/C/ CEIL 3 


00013 

C058 

XP*,CCT,C/ CEU* 


000 14 

CO 5 9 

XP 51 CCT/C/ FLAG 2 


00013 

COftO 

X p 6 > CC T* C / CELLS 


OOui* 

COM 

XpViL'CTiC/ CELL 6 


00017 

006 2 

XP 8 #CCT#C/ CELL7 


OOOjU 

C063 

c ic,ccr,*ooec/ 


OOOZl 

C06* 

C1*>CCTM*OOOC/ 


00022 

C065 

KT*CCt# 0 / 



COf 6 

3*PRGj 30/ 
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C067 

*.C/ Y PEAK BUFFERS 

1 

0003 U 

C06 8 

YPliCCT>C/ FLAd 


00031 

C069 

YP2,CCT,C/ CeLLI 


00032 

C 070 

YP3,L'CT>C/ CELL2 


00033 

C071 

YP*iCCT, 0 / FLAC 2 


0003* 

C072 

YP5 j CC T i C / CELL3 


00033 

C073 

m#C C T,C/ CEU- 4 * 



C07* 

*.CRC/*0/ 


000*0 

C075 

XX1,RCPY,ACC# AC3/ 


000*1 

0076 

XXZ, Jwp,3ASMFG0/ 


000*2 

C077 

FGCjLn,A c r>BA5FfXXl/ 


000*3 

CO70 

jSTjACO/BASEiXX/ 


000** 

C079 

>PULL/AC3/ 


000*3 

CORO 

f PUL/ACO/ 


000*6 

cool 

* R 1 J >0/ 



C0P2 

*.GRG, IOC/ 


ooioc 

CC63 

*,C/ OAT a storage pointer 


CO 5* 

DSPCjncT>*OCC/ 

ooToi 

C0R5 

CSP,CCT,C/ 

TCi 0 . 

oo Jo? 

COPS 

SUN^nCTilOOOflO/ 

Q® 

ooloi 

ooio* 

coe? 

A.C/ GYPC MAIN 


COBB 
C 039 

kx> d ct>o/ 

m,gct, 0 / 

A 


CG90 

I.CRGMOCO/ 



C091 

t.C/ X PEAK INTERRUPT SERVICE 

&§ 

38 

OlAOC 

CQ92 

/ L I j AC 3 #g/ 

OlOCt 

C093 

#RIN# 160/ 

C 1002 

C09* 

jPUSHjAC 0 / 

01003 

C095 

,LG, ACC, BASE, XP1/ 
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0X004 C09& • *0CCH 1/ XPE AK2/ 

0X005 C097 /il/A£C/0/ 1 

01006 0096 /ST,ACC/8ASE/XPl/ 

0X007 0099 /XQ* ACCjBASE/XP?/ 

0X01C CiCO /ST/AC0/8A5E/XP4/ 

OXOH ClOX XPEAk4*lD*AC0,BASE,XP6/ 

0X012 ClC2 #0CC4,2^XP£K4/ 

01011 0103 *CAI#ACO/ 1/ 

01014 CLC4 *SH1>ACQ/4/ 

01015 C1C3 *CA J* AC^/ 1/ 

010X4 C 1C6 /ST,ACC>BA$E*XP3/ 

01017 C iC7 $ J^P*/XP£k'3/ 

0102C C1C8 X?EX4/SHR,AC0,4/ 

01021 C1C9 jsT# ACC, 3flse*XP3/ 

0X022 C110 XPEK5,PULL#ACc/ 

01022 C1U jST,ACCjS*SE/XPZ/ 

01024 C 1 1 2 /PUXL/AC3/ 

01023 CM3 /PbLL/ACO/ 

01024 C114 /RTI,0/ 

01027 C 1 15 XPEAK2,U,AC0,1/ 

01030 C1I6 /ST/ACC/BASE/XPl/ 

01031 C 117 /lG*ACC,BASE/XP3/ 

. 01032 CUB ,ST*ACC,BASE*XP8/ 

f 01033 C119 XP£A*6/U,ACC/8ASE/XP2/ 

oi 01034 C 120 ,BCC*i2/XPEAK7/ 

01033 C121 /CAl/ACO/lZ 

01034 C 122 ,SHR,AC0/4/ 

01037 0123 /C A t* ACO/ 1/ 

01040 C 124 /ST«ACC*9ASEsXP?/ 

01041 C 125 f JMP* # XPEAK8/ 

01042 Cl2ft XPEAK7*5HR/ACC/4/ 

01C43 C 127 /ST/ACC/dASE/XP?/ 

01044 C126 XP£AKS,PCU, ACO/ 

01043 C129 /ST/ACC*BA5E/XP6/ 

01044 C130 /PUU/AC3/ 

OIO 47 CI 3 I /PbLL/ACO/ 

01030 0132 /PTl/O/ 

0133 A,C«Cj12G0/ 

0134 *.0/ V PEAK INTERRUPT SERVICE 
01207 CU5 YPEAK,L.I,AC3/0/ 

01201 C 136 /&1N *169/ 

01202 C 117 / Pb ACC / 

01203 C 1 38 /XC/ACC/BASE/YPl/ 

0X204 CX39 #8CC4,1#YPEAK2/ 

01205 CX40 /XI/ ACC/O/ 

012C6 C 14 X /S t /ACC*BAsC/VP1/ 

0120T C 142 /XC/ACc^BASE/YPS/ 

01? 1C C 1 4 3 ./BCC*j 2/YPEAK3/ 

01211 C 144 /CaI*aC0*1/ 

01212 0145 * SHR ^ ACO* 4 / 

01213 C 146 jCAJ/ACO/4/ 
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01>14 C147 /$U8;AC0;8ASE/YP5/ 

01*13 Cl *8 ;ST> ACC jBASE; YP3/ 

Ol?l* C149 , JKP+,YPEAK4/ 

D1?17 C 1 50 YPEAK3,5h*/AC0>4/ 

0122 ? C 1 5 1 >5L8jAC0iQASE/YP5/ 

01221 C 1 52 iST* ACC,BASE^YP3/ 

01222 C 153 YPEAK4,PLLL/AcO/ 

01223 C 154 ;ST,Aec, 6 ASE;YP 2 / 

01224 C 1 55 ,PIU,AC3/ 

01225 C 1 56 iPLLL^ACO/ 

01226 C 157 ,*11/0/ 

01227 C 158 YPEAK2>U,ACC, 1/ 

0123C C 1 59 #ST,ACC,'3ASr/YPl/ 

01731 C 1 6 0 ;LC,ACC;BA$E;YP2/ 

01732 C 1 6 1 ,BCC+,2,YPEAK5/ 

01 73 3 Clft2 ;CAI,ACC,1 / 

01734 C 163 ;Shft,AC0;4/ 

01733 C 164 jCMjACO,]/ 

01736 C163 -^$Uft , A G0;VP a / , 5U 6, A C cj> > 6 AS £ , Y P2 / 

01737 C 166 /5T#AC0#BA$E/YP6/ 

0124C C 1 6 7 , JKP*iYPEAK6 / 

01741 C 1 6 8 YPEAK5j$hR*ACC»4/ 

01242 C 1 69 ,SlP;ACO;0ASE,YP2/ 

01743 CI 7 O ,ST, ACC, BASF, YP 6 / 

01744 C171 YPEAK 6 ,PLLL,ACQ/ 

01245 C 172 ,ST, ACC, BASF, YP9/ 

01745 C 173 ,PLU,AC3/ 

01747 C 174 ,PUU,AC0/ 

0125C C 175 ,RTt,0/ 

C176 *,CRG,15C0/ 

C177 *.C/ DATA STCRAGE INTERRUPT SERVICE 
0I9CC C 1 78 0S*LC, ACC, BASE, QS P / 

015C1 C 179 ,SH*,ACC,14/ 

01502 C 1 9 0 ,90C*,3,Qs6/ 

01503 Cin ,10, ACC, BASE, XP1/ 

0150* C 1 9 2 ,0OC*,1,OS2/ 

01505“ C i S3 ,IC,ACC, 8 ASE/Xp 6 / 

019C4 C 1 94 >jrr*jLs3/ 

01507 C 185 DS2,LO,ACrj*BA5E,XP2/ 

0131C C 1 86 0S3,IP,AC3,RASE,QSP/ 

015U C 18 7 , ST, ACC, AC 3,0/ 

01512 C 100 , L I , AC C , l / 

0 1 5 1 3 C 189 ; R A DD; AC C; AC 3 / 

02514 C 2 90 ,ID,ACC,BA5E,YP1/ 

01513 C 1 9 1 ,BCC*,1,US4/ 

01916 C 192 ;LC'ACC'BASE/YP3/ 

0ljl7 C 193 ;JFP*;CS 5 / 

0152C C 194 DS4;in,AC0;8ASE#YPZ/ 

01321 C 195 D55,ST,AC0,AC3,0/ 

0152 2 C 1 96 ,LI,ACC;l7 

01323 Cl<?7 ,RADQ,AC0,AC3/ 
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olsiA thi nfc,Acc,SA$E#nT/' ' 

cl^3 iSTiiCC>AC3,0/' v 

01526 C2C0 ,L I,AC2, 1/ 

01327 C2Ci , RAOQ, ACO, AC3/ 

01330 C2C2 ,ST,Ae3,BASE,0SP/ 

01331 C 203 0S6/PUU/AC3/ 

01332 C2C4 ,PULL,ACG/ 

01533 C2C5 ,«Tl,0/ 

C2C6 $.0HG,20C0/ . 

C207 *.C/ GVPC MAIN PROGRAM 

02000 C2C8 *AIN,lO,ACO,BaSE,XX2/ 

02001 C2C9 ,ST,ACC,BASE#XX/ 

02002 C 2 1 0 ,M,AC3,0/ 

02003 C2U ,IC, ACC, BASE, 9SPC/ 

02004 C212 ,ST, aCC,Ba5E,DSP/ 

02003 C 2 1 3 ,Ll,ACC,0/ 

020CA C2U ,$T,ACC,BASEiMX/ 

020C7 C213 ,SFLG,IEN,C/ 

02010 C216 REpLO, ACO, BASE, XM/ 

02011 C217 ,hCC*,l,tfAlN2/ 

02012 C218 ,IC, ACC, BASE, XP8/ 

02013 C219 ,SHL,ACO,A/ 

02014 C220 ,C A I, ACO, 1 / 

02013 C 22 1 ,Acn,AC0,BA$E,XP8/ 

0201fc C222 ,ACD, ACO, BASE, XP6/ 

02017 C223 , ACO, ACO, BASE, XP7/ 

02020 C224 , ACO, ACO, BASE, MX/ 

02021 C225 , ST, ACC, BASE, Mx/ 

02022 C226 ,JMP*;*AIN3/ 

02023 C227 MAIN2,LD,AC0/RASE,XP4/ 

02024 C 2 28 , SHI, ACC, 4/ 

02025 C 2 29 ,CAI,ACO, 1/ 

02026 C230 ,ACO#AcO,BAsE,XP4/ 

02027 C231 ,AC0,AC0iBAS6#XP2/ 

0203C C232 , AC D, AC 0, B A SE, XP3/ 

02031 C233 ,ACD,AC0,8A$E,MX/ 

02032 C234 , ST, ACC, 3 ASE>?1X/ 

02031 C235 MAJN3,8nC*,2,MAlN4/ 
0203* C236 ,C A t , ACO, 1 / 

02033 C237 *ShL,ACO,l/ 

02036 C238 ,rCC*> It Y A INS/ 

02037 C239 , L 0, A CO, 9 AS E , C 14/ 

0204C C240 ,S T, ACO, BASE, MX/ 

02041 C 24 1 , JPP*, MA I N5 / 

02042 C242 MAIN4,SHC, AC0,1/ 

02043 C243 ,RCC*, 2,M A IN5/ 

02044 C 244 ,LC, ACO, BASE, CIO/ 

02045 C245 , ST, ACO, BASE, MX/ 

02044 C 2 46 MAINS, LD, ACC, BASE, YP1/ 
02047 C247 ,BCC*,1,MAIN6/ 

0203C C 24 8 , L I , AC C, 0/ 
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oaosi 

C243I 

.»ado;aco,base,yp3/ 


02032 

C250 

/400>AC0/8ASEi YP6/ 


02033 

C25 1 

4 ST 4 ACC 4 BA 5 E 4 H/ 


02034 

C232 

4 JPP*4FAlM/ 


02033 

C253 

MAIN6/L1/4C0/0/ 


02054 

C 2 54 

,ADC,AC0,BASE # YP2/ 


02057 

C255 

^ACC/AC0^BA5E#YP3/ 


02060 

C256 

aST/ACCa9ASE#H/ 


02061 

C257 

HAlN7un>ACC>DASE*MX/ 


0206 2 

C 2 9 8 

^SL«l,Ac n /8ASE#M/ 


0206 3 

C259 

>ST,ACc/8ASE/M/ 


0206 4 

C 26 0 

^AIM0jPCC4/2,MAIM11/ 


02063 

C26 1 

iCAI/ACO# 1/ 


02O64 

C 262 

,G8,AC0#BA!»e,5lGN/ 


02067 

C 26 3 

^ A I N 1 1 > Ll> Ac3j 0/ 


O2O70 

C264 

* RCUT# 17C/ 


02071 

C265 

jST,ACC#BASEiHT/ 


02072 

C 266 
C 267 

,JHP*,RET/ 

i.C/ PROGRAM SECTION DUMP 


02073 

C 2 68 

DUKPiLD/ACS^aASE^OSPC/ 


02074 

C269 

CP2/8CPY/AC3/ACO/ 


02073 

C270 

/5H*,ACO,14/ 


02076 

C 27 1 

i8CC+/3^0P2/ 


02077 

C272 

/IC» ACOi AC 3 1 0/ 

v> 

1 

021 no 

0273 

4 ShP j AC 0; 1 C / 

oo 

02101 

C274 

4LC*4 AC24PUTC./ 


02102 

C273 

4 JSP4AC2^0/ 


02103 

C276 

4 L C j AC 0/ AC 3/ 0 / 


02104 

C277 

4ANO*4AC04^ASK/ 


0210 3 

C27B 

4 10*4 AC2/ PUTC ( / 


02104 

C279 

4 JSR 4 AC2/ 0/ 


021 07 

C210 

4 U 4 ACC 4 I/ 


02110 

C2P1 

4RA004ACOiAC3/ 


02111 

C 2 82 

4 JNP+4CP2/ 


02112 

CZd3 

PUTC.40CT477131/ 

O'C 

02111 

0234 
C 28 3 
C 206 

^ASK,ucT4377/ 
*.ENC/ 
i • FEE 0/ 
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CIO 

O 0 O 2 C 

C 2 C 44 






CIA 

00021 

C 2 C 37 






□ P 2 

02 Q 7 A 

C 2076 

0 2 i 1 1 





DS 

01500 







0 S 2 

01507 

C 1504 






053 

01510 

C 1506 






OS A 

01520 

C 1515 






053 

01521 

C 1517 






056 

01531 

C i 502 






DSP 

ooioi 

C 150 G 

01510 

01530 

02004 



D$PC 

001 00 

C 2 O 03 

02073 





DLHP 

02073 







FCU 

O 0 OA 2 

C 0 C 41 






H 

00104 

C 2 C 53 

02060 

02062 

02063 



M AIM 

02000 







HA I'.'IO 

02064 







ha mi 

02 <j 67 

C 2 C 64 






HA I M 2 

02023 

C 2011 






HA I M 3 

02033 

C 2 C 22 






HAI *4 

02042 

C 2 C 33 






HA ^5 

02046 

C 2 C 36 

02041 

02043 




HAIV 6 

02035 

C 204 ? 






Ml '.7 

02061 

C 2 C 34 






MAS* 

02113 

C 2104 






HT 

30022 

C 1 524 

02071 





HX ' 

00103 

C 2 CC 6 

02020 

02021 

02031 

02032 

02040 

PLTC. 

02112 

C 2101 

021 C 5 





BET 

n 2010 

C 2 C 72 






S 1 CM 

ooio? 

C 2 C 66 






XPl 

oooio 

C 1 003 

010 C 6 

01030 

C 1503 

02010 


XP 2 

ooon 

C 1 C 23 

01033 

01507 

02027 



XP 3 

00012 

CIC 16 

C 1 0 2 1 

01031 

C 2 O 30 



XPA 

00013 

C 1 C 1 C 

0202 3 

02026 




XP 5 

00014 







XP 6 

00013 

C 1 C 11 

01043 

01305 

02016 



XP 7 

00016 

C 1 C 07 

01040 

01043 

02017 



XPS 

00017 

C 1032 

02012 

02015 




XPEAK 2 

01 027 

C 1 C 04 






XPEAKA 

01011 







XPEAK 6 

0103 ? 





* 


XPEAK 7 

0104 2 

C 1034 






XPEAKS 

01044 

C 1 C 41 






XP 6<4 

0102 J 

C 1012 






XP £><5 

01022 

C 1 C 17 






XX 

00003 

COO 43 

02001 





XXI 

0004 'J 

C 0 CA 2 






XX 2 

000 41 

C 2000 






VPl 

00030 

Cl 20.1 

012 C 5 

01230 

C 1514 

02046 


YP 2 

00031 

C 1223 

01231 

01236 

01242 

01320 

02056 

YP 3 

00032 

C 12 15 

01221 

02057 





02061 
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yneiK2 
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ype*k4 
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00031 
000 3 A 
000 3 3 
012CC 
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0i?i? 
0 1?2 2 
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0 1207 
C 1237 

C1204 
C 1 2 10 
C1216 
C 1 232 
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C298 S 


01214 

01243 
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C290 

t *1? t2* LC/ 


C291 

*.SETC, 6 , 10/ 


C292 

k.JU^A/ 


C 293 

MF,2,LI/ 


C294 

&.SE t C*5,23/ 


C295 

1 • JUH P / R / 


C 296 

*.IF>2,LIN/ 


C2<?7 

l.SETC,5*23/ 


C258 

*,SETC* 6 jC/ 


C 259 

k.S^B,6,4,4/ 


C3C0 

>.JllNP*R/ 


C3C 1 

• 1Fj2«LCI/ 


C 3 C 2 

« SE TC j fr* 1 1 / 


C3C3 

. JIJPP# A/ 


C3C4 

. IF> 2 * LCR/ 


C3C5 

, 5E TO 5> 2300 / 


: C3C6 

•JUPPiK/ 


C 3 C 7 

* IF* 2 * ST/ 


C3ce 

•5E1Cf6« 12/ 


C3C9 

• JU^P> A/ 


C 3 10 

» 1 F# 2>ST 1/ 


C 311 

.SE TC#6> 13/ 


C 312 

•JUPP# A/ 


C313 

•IFj?*ST9/ 


C314 

• SE TC j 3/ 232C/ 


C315 

.JUf'P**/ 

* 

C316 

♦IF#2jRXCH/ 


C317 

.SETC,3,34/ 


C318 

.JU*P,D/ 


C 3 1 9 

• I F* 2# RC PY / 


C3?0 

« SE TCj 5* 35/ 


C321 

. Ju^n/ 


C 322 

, IF,2,PLSH/ 

o P 

C323 

.S€TC,4,20/ 

td 

C324 

• jUt'PfC/ 

*p o 

C 32 3 

. IF>2>PtLL/ 


C 32 6 

fc - SETC# 4^ 21/ 

§ > 

C327 

,JUKP,C/ 


C 328 1 

,IF,2,XCHRS/ 


C329 1 

•SETCt 4j 23/ 


C330 1 

.ju*p,e/ 

fe c 

C331 1 

.!Fi2>PtSHF/ 

S fe: 

C332 1 

.SEK«3/1/ 

Sir 

C 3 3 3 1 

•JUFPaL/ 


C334 < 

• IF#2aPLUF/ 


C335 1 

■SETC/3/5/ 


C 3 36 1 

.jurpu/ 


0337 | 

, 1 F, 2j ADD/ 


0338 1 

• SETC#.6# 14/ 
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C3<*0 ».lF,2ip*DD/ 

C 3 A 1 *.$£TC,5,22AC/ 
C3A2 *.J ! JKP;K/ 

C3A3 |.IF,2,CAI/ 

C3AA »,5ETC,5,2*/ 

C3A5 S.JUh'PjB/ 

C3A6 *.IF,2/RA0D/ 

C 3A 7 4.5ETC>5,30/ 

C3A8 4,*IU*P,P/ 

C3A? ».IF,2,SLB/ 

C350 4 , SETCj 6* 13/ 

C 35 1 * . JUf'Pj. A/ 

C 35 2 4.lF>2,nSU9/ 

C353 *.sETC/5/Z2ftO/ 
C35A 4, JUMP, K/ 

C 3 5 5 I* lF/2jMf‘Y/ 

0356 I « SETCj 5*22^0/ 
C357 *.JU»'P,K/ 

C 35 8 4,IF«2#D]V/ 

C 359 *.SETCi3#2220/ 
C360 4 « JUKP/ k/ 

C 3 1 1 *'IF«2/A.\0/ 

C362 4,SETC,6/H/ 

C363 I.JUf'Ptf,/ 

C3ft* 4.IF,2,M*JD/ 

C3<55 4,SETC#5/ 37/ 

C366 *,JU^P>0/ 

C367 $ « IF# 2#OR / 

:aea 4 .sE7cj6>i5/ 

C3F9 I.JU/PiG/ 

C 370 4.IF,2,RxOR/ 

C37i 4 » 5 F TC j 5* 3ft/ 

1372 i.JO^PjP/ 

■2 373 4.iF,2,SPT0 I T,3/ 
:37A 4 , SETC* A/ 162/ 
,0375 I.Ju^P>J/ 

2376 f . IF> 2^CL»BIT>4/ 
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Direct-axis and cross -axis digital transfer functions 

Mechanized direct-axis and cross-axis digital transfer 
functions 

General analog transfer function 
General digital transfer function 
Dummy analog transfer function 
Iteration index 
Pickoff gain factor 

Compensation gain factor 

Torquer gain factor 

Total digital torque 

Components of computer digital torque resulting from 
X- and Y- pickoff signals 

Analog rebalance torques 

Gyro spin speed 
P recessional torques 

Laplace Operator 
Total torques 

Digitized gyro pickoff angles 
Sampled X-axis pickoff angle 

Sampled Y-axis pickoff angle 
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Delay operator 

Digitized gyro pickoff angles 

Gyro pickoff angles 

Iteration period 
Angular rate inputs 

Analog filter break frequencies 

Dummy analog filter break frequencie 



